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Abstract 



Phase transitions of small isolated systems are signaled by the shape of the 
caloric equation of state e* (T) , the relationship between the excitation energy 
per nucleon e* and temperature. In this work we compare the experimentally 
deduced e*(T) to the theoretical predictions. The experimentally accessible 
temperature was extracted from evaporation spectra from incomplete fusion 
reactions leading to residue nuclei. The experimental e*(T) dependence ex- 
hibits the characteristic S-shape at e* = 2 — 3 MeV/A. Such behavior is 
expected for a finite system at a phase transition. The observed dependence 
agrees with predictions of the MMMC-model, which simulates the total ac- 
cessible phase-space of fragmentation. 
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In the macroscopic physics phase transitions are usually defined by a divergence at the 
critical temperature, for example in heat capacity c = de*/dT thd , where e* and T thd is 
the excitation energy and the thermodynamic temperature. This is corresponding to the 
well known finding that at a first order phase transition temperature stays constant while 
additional energy is pumped into the system. This picture becomes different if we deal with 
finite and isolated systems like nuclei. Due to conservation of mass, charge and especially 
total energy the signal of a first order phase transition is given by an "S-shape" in e*(T), 
called the caloric equation of state (CES), as shown in the figure [j] for a decaying nucleus. 
Pictorially speaking we find that the system is cooling down with rising excitation energy at a 
first order phase transition in a finite system [l|-|5j . For a finite and isolated (microcanonical) 
system the heat capacity is no longer a positive definite quantity. At first-order phase 
transitions it has two divergences (instead of one for the infinite matter). In the region 
between the two poles it becomes a multi-valued function. 

The signal of fig. [I] can be obtained only in a microcanonical description which takes into 
account the strict mass, charge and energy conservation. This behavior of a fragmenting 
nuclear system at a phase transition is due to the opening of new decay channels, i.e. 
the population of additional regions of phase space Q(e*) (e* is here the excitation energy 
per nucleon) 0|§. In case of fig. [I] it is connected to the onset of IMF (intermediate 
mass fragment) emission and the new phase space associated to IMF. The specific entropy 
s(e*) = ZwS jl 6 \ N the number of nuclei, rises then higher compared to a normal Fermi-gas. 
It is this strong rise of entropy that leads to an anomaly in the CES by the relation 

1 -«0 = 5gi . (i) 



A recent experimental observation |n|,fnj showing a structure in e*(T)-curve fueled the 
discussion about the appearance and the measurement of a phase transition. It is known that 
apparent temperature is sensitively dependent on the mass of the source M. We suppose 



together with |TB| , |H| and in opposite to jnj that the curve shown in ref. [IT| is just the 
effect of changing mass of the source without undergoing any phase transition. Another 



discussion can be seen in [HJ. Here we concentrate on a different set of experimental data 
deduced from ref. (IT! 



We are going to perform a comparison of an experimentally obtained "S-shape" in e*(T) 
with theoretical predictions of the Berlin - microcanonical statistical multifragmentation 
model MMMC which simulates the phase space Q(e*) for decaying nuclei. We are going to 
describe the model in some detail at a later stage. MMMC predicts two phase transitions 
in nuclear fragmentation |J. The phase transition at lower excitation energy at e* ~ 2 — 3 
MeV per nucleon was shown in figure [l]. A second phase transition at higher e*, which is 
not the subject of this treatise and not shown here, is due to the true multifragmentation. 

Another similar statistical fragmentation model, the Copenhagen model SMM also 
predicts phase transitions. Since SMM has some mixed microcanonical-canonical features 
and has a varying freeze-out volume || it produces a slightly different signal of a phase 
transition compared to MMMC. 

The thermodynamical temperature T thd , equation ([[]), cannot be accessed directly in 
an experiment. For the experimental comparison we need to find a related quantity which 
would keep the information on the behavior of T t hd [PIP1- Such a quantity, which we call 
apparent temperature T app is thus not a temperature in the sense of thermodynamics. 

In this work we show experimentally accessible "S-shapes" of the CES e*(T app ) extracted 
from incomplete fusion reactions resulting from 701 MeV 28 Si + 100 Mo |17| . We plot e* vs. 
T app , where T app (the apparent temperature) is the slope of the raw evaporation spectra. 

The details of the experiment and the extraction of the needed parameters can be found 
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17| , |18| . Here we outline some of the important features. Heavy evaporation residues 
were detected at forward angles, therefore this experiment does not probe multi-fragment 
final states. Charged particles (including IMFs) and neutrons were detected in concentric 
47r-detectors. The excitation energy of the source was deduced from linear momentum 
reconstruction. The raw spectra of protons, deuterons, tritons and alpha particles were fitted 
with a three moving source prescription. The data at backward angles are well described by 
a surface-evaporating Maxwellian moving source: 



where E^m is the center of mass kinetic energy of the particles, B the Coulomb barrier and 
T app , which is the slope of the raw spectra, is the desired apparent temperature. 

Figure ^| presents the excitation energy per nucleon e* versus T app for protons, deuterons, 
tritons and alpha particles. This representation of the data exhibit two noteworthy trends. 
The first trend concerns the general shape of these curves and the second is the horizontal 
displacement (along the T app axis) as one progresses from protons to deuterons to tritons 
and alpha particles. We shall focus on the first observation although the second observation 
is also of interest and we shall briefly discuss it also. We find that all the four curves for 
different particles show an " S-shape" in the expected region of excitation energies (compare 
to fig. P, but no backbending is seen. Here one needs to keep in mind that the experimental 
data points correspond to sources which are slightly changing with the excitation energy. 
From lowest to highest energy the mass is growing from 105 to 122 nuclei and the charge 
from 47 to 54. We expect this change in mass and charge to smear out the " S-shape" with 
a backbending shown in figure [l|. 

Next we are going to perform a comparison to the MMMC-model simulation. 

The MMMC-model assumes that the compound system fragments quite early but the 
fragments remain stochastically coupled as long as they are in close contact. Consequently, 
the system is equilibrated inside a freeze-out volume. The size of this volume, which is 
a simulation parameter of MMMC is in our energy region at about 6 times the normal 
nuclear volume. This corresponds to an average maximum distance between the fragments 
of ~ 2fm. Here the nuclear interaction between the fragments drops to the point that 
subsequent mass exchange is unlikely. Then the fragments (which can be in excited states) 
leave this volume and may de-excite as they trace out Coulomb-trajectories. The ensuing 
formation of fragments is determined by the accessible phase space which is sampled with 
the Monte Carlo method using the Metropolis importance sampling. 

The experimental analysis of the data provides the values of the mass A, charge Z, 



excitation energy E*, and angular momentum L of the source fll7[ , which are the input 
into the MMMC simulation. The only simulation parameter of the model, the freeze-out 
radius Rf was taken as its standard value of 2.2A 1 / 3 fm, this means that we simulate a 
phase transition at constant volume. The results of MMMC calculations, performed with 
these input values, were subjected to a software filter of the experimental set-up which, most 
importantly, selects only those events with one big residue. The mass of the residue was 
chosen to be A res > 90, which is close to A res estimated from the experimental data (the 
experiment did not directly measure the mass of the residue). 

Figure ^| shows a comparison of the e*(T app ) curves extracted from the experimental data 
for protons and a-particles to the e*(T app ) dependence deduced from the MMMC-mode\ 
using its standard parameters. Also the experimental uncertainties for the proton and 
alpha curves are given. The horizontal bars give the statistical uncertainty to extract the 
temperature (slope) from the experimental raw spectra. The vertical bars (here only given 
at the lowest and highest proton or alpha point) indicate the systematic difference of the 
excitation energy extracted by the "top-down" resp. the "bottom-up" procedures employed 
in ref. |T7J . The two alternative methods lead essentially to an up or down shift of the CES 



curve without changing the main structure of the curves. The similarity of the shapes of 
the experimental and simulated CES e*(T app ) is quite evident. The differences between the 
shapes of these curves and the parabolic dependence (dotted curve) expected for a simple 
Fermi gas is clearly seen indicating that some additional degrees of freedom, which are 
apparently included in the MMMC-model, become significant in this energy range. 

The theoretical value of T app was extracted from fitting, as was done for the raw ex- 
perimental spectra. It is worth noting that calculated temperatures T app , extracted from 
the Maxwellian fits, are close though not identical to the unique thermodynamical temper- 
atures T t hd from the equation f[[ as can be seen from comparing figure [I] and figure |3]. The 
curve in figure |1] is calculated for the mass and charge corresponding to the highest value of 
experimental energies, but for the whole energy range. 

The values of Rf and A res do not influence the general shape of the calculated e*(T app ) 



curves. However, the e*(T app ) curves shift along the T app - axis if different values of these 
parameters are used. The shifts produced by reasonable changes in A res are larger than 
those produced by reasonable changes in Rf. We checked that the anomaly in the CES is 
not due to the changes of the angular momentum from L = 18.2 to 48.8ft. It exists also at 
L = 0. 

While the similarity of the shapes of the experimental and simulated CES e*(T app ) is 
good for p's and a-particles, significant differences exist. The simulated curves for deuterons 
and tritons (not shown here) have the same shape but are shifted towards lower values of 
T app . The higher T app values of the experimental deuteron and triton spectra might be an 
indication that the production of these less bound fragments might occur in an earlier hotter 
stage of the reaction, an ingredient not included in the MMMC . 

We have also compared the multiplicities of neutrons M n , protons M p , deuterons M d , 
tritons M t and alpha M a calculated with MMMC with the experimental values. The total 
number of the evaporated particles is the same in the calculation and experiment. The 
model overestimates M p by approximately a factor 1.3 and underestimates M a by a factor 
of 2. The MMMC calculation reproduces M n and M t at all values of e*. On the other hand, 
the values of M d are not reproduced by the model calculation, which systematically predicts 
values which are too high. In this context one should keep in mind: In MMMC we treat 
deuterons as spherical nuclei with normal nuclear matter density. This may overestimate 
their stability. 

The experimental data also suggest an association between the onset of significant IMF 
production and the upswing in the e*(T app ) dependence. This is seen in figure |5] where the 
measured absolute IMF-multiplicities (Mimf) associated with the experimentally selected 
events which produce a large residue are compared to the absolute M IMF of the correspond- 
ing events from the model calculation. (The smooth dependence of Mimf on the excitation 
energy underlines the high statistical quality of both experimental and simulated data.) 
Both the data and the calculations exhibit a dramatic increase in Mimf for e* between 
e* ~ 1.5 and ~ 3MeV/A. The primary difference between the experimental data and the 
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model predictions is that the M/mf's rise at higher e*- values for the model predictions than 
they do for the experimental data. Another point worth noting is that the values of Mjmf 
in these decay-channels are less than 0.1. In other words, most (more than 90%) of the 
selected events (both in the experiment and in the model calculation) have no IMF. 

Besides the overall agreement in the experimental an theoretical CES some uncertainty 
of the interpretation remains. There are too few experimental data points outside the "S- 
shape" region to interpret the data unambiguously. Further there is a small change in the 
charge and mass of the source for different excitation energies. Therefore it is desirable to 
perform a similar experimental analysis covering a larger range of excitation energies and 
selecting a strictly constant A source and Z source for the whole energy range. 

Finally, we would like to point out that the extracted temperatures T app for p and a- 
particles in related earlier experimental work by the Texas A&M-group [pi]] , exhibit trends 
similar to those presented here. This data is for similar masses of compound nucleus ranging 
from 109 to 128 and in the same energy region. We plot the data for the apparent tempera- 
tures T app (e*) for alphas in figure f|. The proton data (not plotted here) show also a similar 
backbending. Despite the large error-bars this shape anomaly had already been noted and 
was well reproduced by MMMC in fig. 13 of ref. M. 



In this paper we have shown that a strong anomaly exists in the shape of the experimental 
CES e*(T app ) for the apparent temperature. This "S-shape" suggests that the relevant phase 
space becomes enlarged in the region of e* = 2 — 3MeV/A. In terms of thermodynamics this 
is associated with a phase transition for this isolated, strongly interacting quantal system. 
MMMC-model reproduces the general shape of the experimental e*(T app ) curve at right 
excitation energies. This supports the hypothesis of strong stochastic mixing of the various 
fragmentation channels and the statistical equilibration at freeze-out. As the production 
of intermediate mass fragments increases dramatically in this region of excitation energy 
we associate the "S-shape" to the additional phase space opened by IMF production. This 
association is also supported by the results of the MMMC calculation where the " S-shape" 
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in e*(T) is seen in the energy region of strongly increasing Mjmf- 

More to the point the "S-shape" in the caloric equation of state associated with IMF 
production is even seen in the evaporation spectra of different particles in events which in 
more than 90% of the cases have no IMF. In addition, the calculations are rather insensitive 
to variations of its basic model parameter, the freeze-out volume within broad limits || and 
thus no adjustment of this parameter was necessary to reproduce the general shape of the 
experimental caloric equation of state. Prior to this work the primary evidence of the validity 
of the concept of stochastic coupling of two moving nuclei in proximity was the finding of a 



strong pre-barrier surface friction in deep inelastic collisions ||19|| . An experimental support 
for the equilibration hypothesis was given also at higher excitation energies in |20| . 



While this work adds weight to the argument that strong stochastic mixing and equi- 
libration exists up to rather extended configurations of the fragmented source, we do not 
believe that the issue is closed. For while the general shape of the caloric equation of state 
was reproduced by the MMMC model, differences exist with particle type which may imply 
the existence of a dynamics or a (mean) sequence, which is not dealt with by the single 
freeze-out configuration of the MMMC model. 

Taking all the findings together, the anomaly in all four spectra (proton, deuteron, triton, 
and alpha) at the same excitation energy as predicted by MMMC and also the earlier data 
by the Texas A&M-group (fig. |2|, |3] and |j), we see a strong support for the significance 
of our interpretation of the "S-shape" in the e*(T app ) cLS cL SI gnal of a phase transition 
in nuclear fragmentation. The transition is from pure evaporation to asymmetric fission, 
which is associated to the onset of IMF emission. Nevertheless, additional experimental and 
theoretical confidence is desirable. 

O.S. is grateful to GANIL for the friendly atmosphere during her stays there. This work 
is supported in part by IN2P3/CNRS. 
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FIGURES 
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FIG. 1. Theoretical first order phase transition predicted with MMMC. e*{T t hd) is calculated 
for A source = 122, Z source = 54 and zero angular momentum. The computational error-bars are of 
the size of the symbol. 
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figure 2 

FIG. 2. Experimental excitation energy per nucleon e* versus apparent temperature T c 
p, d, t and a. 

The error-bars are given in fig.[|. 
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figure 3 



FIG. 3. Experimental and theoretical (with MMMC) caloric equation of state, e*(T app ) for 
p and a. The horizontal error-bars give the statistical uncertainty to extract the slope from the 



raw spectra in ref. 17]. Different methods to determine the excitation energy lead essentially to a 
parallel up or down shift of the curves by the amount indicated by the vertical bars at the lowest 
and highest data point. The dotted curve is a Fermi-gas calculation (E = T 2 /10). 
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FIG. 4. The caloric equation of state for a using data from ref. [21]. 
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FIG. 5. IMF multiplicities per event from the experiment and the MMMC-model 
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